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Landau-Zener tunneling of a single Tb^+ magnetic moment allowing the electronic 

read-out of a nuclear spin 
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(Dated: April 25, 2013) 

A multi-terminal device based on a carbon nanotube quantum dot was used at very low tem- 
perature to probe a single electronic and nuclear spin embedded in a bis-phthalocyanin Terbium 
(III) complex (TbPc2). A spin- valve signature with large conductance jumps was found when two 
molecules were strongly coupled to the nanotube. The application of a transverse field separated 
the magnetic signal of both molecules and enabled single-shot read-out of the Terbium nuclear spin. 
The Landau-Zener (LZ) quantum tunneling probability was studied as a function of field sweep rate, 
establishing a good agreement with the LZ equation and yielding the tunnel splitting A. It was 
found that A increased linearly as a function of the transverse field. These studies are an essential 
prerequisite for the coherent manipulation of a single nuclear spin in TbPc2 . 

PACS numbers: 75.50.Xx, 81.07.Nb, 81.07.Ta, 85.75.-d 
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I. INTRODUCTION 

The detection and manipulation of single localized 
spins is a rapidly developing field of nanoscience, en- 
abling the creation of quantum bits and memories. 
Among the large variety of spin systems, electron spins 
confined in quantum dots^^, impurities,^ defects in 
semiconductors,^ and nitrogen- vacancy center^^^^^ are 
promising candidates for spin-based quantum computa- 
tion. However, control of the decoherence induced by 
the environment^^ and control of the coupling between 
different spins are still experimental challenges. Single- 
molecule magnets (SMMs) offer an interesting alternative 
for the embedding of spins in a controllable and repro- 
ducible environment. ^^^ 

A SMM consists of an inner magnetic core, contain- 
ing in general several magnetic atoms, and an outer 
non-magnetic shell, whi ch is chemically tailored through 
the addition of ligands.l^I^ At low temperatures, the 
magnetic core of an SMM behaves like a single mag- 
netic moment, which can attain spin ground states up 
to S = 83/2.^^ SMMs are characterized by a broad range 
of quantum effects, ranging from quantum tunneling of 
magnetization (QTM)^^ ^^ to Berry-phase interference.^^ 
Their spin coherence time is expected to reach the mi- 
crosecond range^^^^, and the perspective of entangle- 
ment between SMMs^ makes them good candidates 
for quantum spintronics^^ and quantum computations.^^ 
Among the large variety of known SMMs, rare earth 
based SMMs are of particular interest because the high 
intrinsic anisotropy of rare earth elements can lead to 
larger zero-field splittings than for transition metal based 
SMMs.^^ Furthermore, the strong hyperfine interaction 
of rare-earth elements leads to well separate resonant tun- 
nel transitions for each nuclear spin state. ^^ Because of 
the excellent intrinsic isolation of nuclear spins, they are 
promising candidates for spin quantum bits.l^^^ 



In this paper, we demonstrate the probing of single 
electronic and nuclear spins in rare earth based SMMs 
adsorbed onto the sidewall of a carbon nanotube. The 
present study focuses on the so-called Terbium "Dou- 
ble Decker" or bis-phthalocyanin Terbium (III) complex, 
referred to as TbPc2 in the following, which was synthe- 
sized using the method reported in Ref.^. 



II. SPIN HAMILTONIAN AND 
LANDAU-ZENER TUNNELING IN TbPcs 

The TbPc2 complex is composed of two organic ph- 
talocyanins (PC2) coordinating a single Terbium ion in a 
square antiprismatic geometry [Fig. 1(b)]. The 4f^ elec- 
tronic shell of the Tb^+ leads to a total magnetic mo- 
ment 7 = ^ + 5 = 3 + 3. In addition, the Tb ion has 
a nuclear spin / = 3/2. The 7r-type electronic struc- 
ture of Pc2 permits coupling between the Tb and the 
environment j^^lMl j^ order to calculate the Zeeman dia- 
gram of TbPc2 [Fig. 1(a)], we used the ligand field Hamil- 
tonian expressed in the Stevens operator formalism)^ 

Uif = aA%Ol + ^{AlOl + AlOt) + ^A%Ol (1) 



and the Zeeman Hamiltonian: 



Uz 



eeman 



5jMoMbJ-H, 



(2) 



where the O^ are the Stevens operators^, A^ the lig- 
and field parameters from Ref.^^, a, P and 7 the Stevens 
constant tabulated in Ref.™ gj the gyromagnetic fac- 
tor of Tb, J the electronic magnetic moment operator of 
Tb, and H the applied magnetic field. Numerical diago- 
nalization of Hzeeman + y^if gives | J = 6, J^ = ±6) as a 
ground state doublet, isolated from the first excited state 
I J = 6, J^ = ±5), by about 600 i^ [Fig. 1(a)]. Therefore, 
the Tb magnetic moment is aligned with the quantiza- 
tion axis (easy axis) within the low temperature and low 
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FIG. 1. (Color online) (a) Zeeman diagram of the ground 
state J = 6. I J = 6, J^ = ±6) is the low lying ground state 
doublet, which is separated by about 600 K from the first 
excited state \J = 6, J^ = ±5). The two mechanisms for 
moment reversal are quantum tunneling of the magnetization 
(QTM) occurring at zero magnetic field and the direct tran- 
sition (DT) occurring at finite magnetic field, (b) Structure 
of the TbPc2 complex. The black, blue, white and purple 
atoms correspond to C, N, H, and Tb, respectively . The 
magnetic easy axis and the hard plane are respectively or- 
thogonal and parallel to the phtalocyanins. (c) Zoom on the 
avoided level crossing label by QTM in (a). A is the tunnel 
splitting, and Pqtm is the tunnel probability, (d) Zeeman 
diagram of the ground doublet | Jz = ±6) split by the strong 
hyperfine interaction of Tb. The color code illustrates the 
different nuclear spin states and the circled intersections are 
avoided level crossings. 



magnetic field regime (T< IK^ jiqH < 2 T) and is consid- 
ered as an Ising-type spin in the following. The presence 
of a transverse anisotropy in the ligand field Hamiltonian 
(A4O4) couples the two states \Jz = ±6), giving rise to an 
avoided level crossing with a tunnel splitting A ?^ 1 fiK 
[Fig. 1(c)]. As a result, quantum tunneling of the mag- 
netization (QTM) can occur in TbPc2 when the mag- 
netic field is swept through this avoided level crossing. 
The tunneling probability is given by the Landau- Zener 
equation?^ 



Pqtm = I - exp 



-^A^ 



^hgjiiB \Jz\l^Q 



dt 



(3) 



where dH/dt is the magnetic field sweep rate and J^ = 
±6. When the electronic moment does not tunnel at the 
avoided level crossing, it ends up in the excited state and 
can relax to the ground state via a direct transition (DT). 
This process involves a spin-phonon interaction in which 
the Zeeman energy is released via a phonon. It occurs 
in the presence of a magnetic field typically larger than 
Q2T.I301 

An important property of the Tb ion is the presence 
of a / = 3/2 nuclear spin. The strong hyperfine interac- 
tion with the electronic moment splits each ground state 
\J = Q^Jz = i6) into four substates. We considered 
the following spin Hamiltonian containing the ligand field 
interaction, the Zeeman term, the hyperfine interaction 
and the quadrupolar term: 



n = Ulf + Uze 



AI.J + P(4^-i(I + l)I) (4) 



where A is the hyperfine constant, P the quadrupolar 
constant and I the Tb nuclear magnetic moment opera- 
tor. The corresponding Zeeman diagram of the ground 
multiplet is plotted in Fig. 1(d). Among the sixteen level 
intersections, only the four encircled intersections cor- 
respond to avoided level crossings. Tunneling can take 
place at one of these avoided level crossings. It changes 
the electronic magnetic moment by A J^ = ±12 but con- 
serves the nuclear spin state. As a consequence, the mea- 
surement of the field position of QTM in TbPc2 enables 
a direct read-out of the Tb nuclear spin state. 



III. SAMPLE FABRICATION AND COULOMB 
BLOCKADE 

QTM was studied in individual TbPc2 complexes, lat- 
erally coupled to a carbon nanotube quantum dot. Cat- 
alyst islands of nanoporous alumina with Fe/Mo were 
first patterned on a Si02 surface with a metallic back- 
gate by optical lithography using L0R3A resist. The 
nanotubes were then grown in a Firstnano CVD oven at 
750° C, whereby methane was used as a carbon source. 
Source, drain, and side gate electrodes were patterned 
by aligned electron-beam lithography, defining 200 nm 
long CNT junctions. The residual resist was removed 
by annealing the sample at 300° C under Ar flow. The 
bis-phtalocyanin Tb(III) complexes were functionalized 
with one pyrene and six hexyl groups in ord er to im- 
prove their grafting onto the carbon nanotub d^^ * ^^ * and 
prevent their crystallization. They were suspended in 
dichloromethane and drop-casted onto the nanotubes. 
Samples with large source-drain resistance (>100 kl] at 
room temperature) were micro-bonded and measured in 
a dilution fridge with a base temperature of about 40 
mK. The electronic temperature of the source and drain 
electrodes was estimated to be about 150 mK. The setup 
was equipped with two magnets, generating up to 1 T 
in two orthogonal directions. Magnetic moment rever- 
sal of the TbPc2 was detected by means of conductance 
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FIG. 2. (Color online) Stability diagram of the device shown 
in the inset. The red lines are degeneracy points in the 
Coulomb blockade regime. Coupling between the leads and 
the quantum dot can be tuned by moving along one of these 
lines, by means of the back-gate and side-gate voltages (Vbg 
and l^g, respectively). Inset: Scanning electron micrograph 
of the studied device. The green electrodes correspond to the 
source and drain, and the blue electrode corresponds to the 
side gate. 



measurements through the carbon nanotube, as proposed 
previously. ^^ The differential conductance was measured 
as a function of side-gate and back-gate with an Adwin 
real-time acquisition system, programmed in a lock-in 
mode (100 /iV and 67 Hz) at zero bias. Fig. 2 presents 
the differential conductance measurements of this four- 
terminals device, shown in the inset of Fig. 2. The con- 
ductance lines correspond to the Coulomb peaks (charge 
degeneracy state), and the zero-conductance region cor- 
responds to Coulomb blockade. Both the side-gate and 
the back-gate are used to tune the chemical potential of 
the quantum dot, and to modify its coupling with the 
leads and the TbPc2 molecules. 



are needed to fully understand this effect. 

For the present study, the conductance of the device 
in Fig. 2 was recorded as a function of the magnetic field 
for each non-zero conductance region. The spin-valve 
feature presented in Fig. 3(a) was observed at Vbg = V 
and Vsg = V. Along the corresponding degeneracy line, 
the signal vanished. One reason for this gate dependence 
could be that the coupling between the molecule and the 
nanotube, mediated by the tt type electrons of the Pc2, 
was changed by the electric field. Indeed, it was shown by 
Lodi Rizzini et al^ that the exchange coupling between 
the Tb and the environment is very sensitive to the oxi- 
dation state of the Pc2. This possible explanation has to 
be strengthened by further experiments and theoretical 
calculations. 

This spin dependent effect is best visible when the 
two molecules have different coercive fields: in this case, 
each magnetic moment reversal induces a sharp change 
in the conductance of the device at a different field. In 
the present sample, one of these molecules (referred to 
as molecule A in the following) had a QTM probabil- 
ity close to one, inside the experimentally resolved sweep 
rate window, while the other molecule (molecule B) had 
a sweep rate dependent probability. It is not surpris- 
ing that both molecules had a different tunnel splitting 
amplitude, since any deformation of the molecule (e.g. 
a slight twist between both Pc) generates a significant 
change in the longitudinal and transverse anisotropy.l^ 

For instance, for rather fast field sweep rates, A ex- 
perienced QTM while B experienced a DT. The corre- 
sponding magneto-transport measurement are presented 
in Fig. 3(a). For slow field sweep rates, the mag- 
netic moment of both molecules reversed by tunneling 
[Fig. 3(b)] and it was therefore difficult to distinguish be- 
tween molecule A and B. However, because the molecules 
had their easy axes aligned in two different directions 
[Fig. 4(a)], with a constant transverse field Hj_, orthogo- 
nal to the easy axis of molecule B, we could discriminate 
between the QTM positions of molecules A and B. In- 
deed, as depict in Fig. 4a, the longitudinal field felt by 
molecule A is H±SinO + H\\CosO while the field felt by 
molecule B is simply H\\^ where is the angle between 
the easy axes of molecules A and B. More details con- 
cerning the magnetic field alignment are presented in the 
following section. 



IV. THE SPIN- VALVE FEATURES 

We recently discovered a surprisingly strong interac- 
tion between the charge carriers of a carbon nanotube 
and TbPc2 molecules grafted onto its sidewall.1^ We 
found that this interaction gives rise to a spin-valve sig- 
nature in the magneto- transport measurements: when 
the magnetic moment of two molecules are aligned in the 
parallel configuration the conductance is maximum while 
in the antiparallel configuration the conductance is mini- 
mum. Even so two theoretical studies proposed an inter- 
pretation of the underlying mechanism, ^^ ^^ more studies 



V. DETERMINATION OF THE EASY AXIS OF 
MAGNETIZATION 

In order to determine the easy axis of magnetization of 
molecule B, the spin- valve effect was measured as a func- 
tion of the applied magnetic field angle. For each angle 
we measured the magnetoconductance for the trace (from 
-IT to +1T) and the retrace (from +1 T to -IT). A 
relative fast sweep rate of 50mT.s~^ was used to have a 
high probability for a DT of molecule B. Fig. 4(b) shows 
the difference between retrace and trace for positive mag- 
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FIG. 3. (Color online) (a) Magneto-conductance recorded at 
a sweep rate of 50mT.s~^. The blue curve corresponds to 
the measurements made from — 0.35T up to +0.35 T, and 
the red curve corresponds to the ones made from +0.35 T 
down to —0.35 T. Each of the sharp changes in conductance 
corresponds to the signature of a single TbPc2. The reversal, 
which takes place under a low magnetic field, is due to QTM, 
whereas that which takes place under a higher magnetic field 
is due to DT. (b) Magneto-conductance recorded at a sweep 
rate of 20mT.s~^, the colour code is the same as for (a). 
The sharp conductance changes occur close to zero magnetic 
field and are the consequence of both molecules experiencing 
QTM. 



netic fields and the difference between trace and retrace 
for negative fields. The red areas between —70° and +70° 
correspond to curves similar to Fig. 3(a), whereas the re- 
gion between +70° and +110° correspond to a situation 
where the magnetic field is not strong enough to allow a 
DT [see Fig. 3(b) of Ref.^^J. From this characterization, 
the easy axis of molecule B was found to be orthogonal 
to the nanotube axis As a consequence, the transverse 
field was therefore applied along the nanotube axis. 



VI. ELECTRONIC READ-OUT OF A SINGLE 
NUCLEAR SPIN 

Figure 5a shows measurements recorded with a con- 
stant transverse field of 0.35 T and presenting only QTM 
features, the ones presenting DT being rejected. The 
conductance suddenly decreases, between the magnetic 
field values of -50 mT and +50 mT (QTM of B), and 
then increases in all cases above approximately +50 mT 
(QTM of A). This measurement was repeated 3500 times 
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FIG. 4. (Color online) (a) Schematic of the relative orienta- 
tion of the two molecules on the nanotube. We define the lon- 
gitudinal if II and transverse H^ with respect to the molecule 
B which owns its easy axis perpendicular to the nanotube axis 
(see Fig. 4(a)), and Q the angle between the easy axis of B 
and A. (b) Anisotropy of the direct transition. The difference 
between the trace (from — 1 T to +1 T) and the retrace (from 
+ 1T to —IT) is recorded as a function of angular orienta- 
tion. The locus of the direct transition reversal is indicated by 
dashed lines. It is important to note that the reversal close to 
the 90° orientation cannot occur because of the field limita- 
tions imposed by the magnetic field coils in the experimental 
setup. 



with a lOOmT.s"^ sweep rate. Whenever DT failed to 
occur, the longitudinal position of the QTM in molecule 
B was stored in the histogram plotted in Fig. 5b: four 
peaks emerge with a FWHM of approximately 10 mT, 
and a mean peak-to-peak separation of 25 mT. 

In order to explain these results, we compare the posi- 
tion of these peaks with the Zeeman diagram. Fig. 5(b) 
shows the very good correspondence between the four 
peaks and the avoided-level crossing of the Zeeman di- 
agram. This diagram is slightly different from the one 
presented in Fig. 1(c) since we took into account that 
the easy axis is not lying exactly in the plane {H\\,H_\_). 
As evident from the comparison between the histogram 
and the diagram, each of these peaks can be attributed 
to a particular nuclear spin state. 

Similar results have recently been demonstrated at the 
single molecule level by Vincent et al.^, in a molecular 
transistor configuration. It is important to note that in 
the present case the FWHM is larger than in the case 
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FIG. 5. (Color online) (a) Magneto-conductance curves 
recorded under a 0.35 T transverse field, at a lOOmT.s"^ 
sweep rate. The relative QTM position for molecule A and 
B are clearly split, the angle between their easy axes is then 
estimated around 15°. (b) Histogram of the QTM position 
of molecule B, for 3500 consecutive traces recorded under the 
same conditions as those applied in (a). The position of those 
peaks corresponds to the avoided level- crossing shown on the 
Zeeman diagram of the ground doublet Jz = ±6 split by the 
hyperfine coupling with the nuclear spin / = 3/2. 



of Vincent et al. leading to a lower fidelity in the single- 
shot read-out measurement. One reason for this could be 
that the current tunneling through the carbone nanotube 
is interacting more strongly with the TbPc2 molecules. 
Nevertheless, one advantage of the present device is the 
very large variation of conductance induced by the spin 
reversal (200% in the present case versus 1% in the work 
of Vincent et al.^^, which makes the measurement very 
easy since it does not require any specific filtering (phys- 
ical or numerical). 

In order to confirm that only one QTM position ex- 
ists per nuclear spin state, we measured the tunneling 
probability as a function of the sweep rate and com- 
pared it with the Landau- Zener theory. 100 magneto- 
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FIG. 6. (Color online) QTM probability as a function of the 
inverse of the sweep rate, for different transverse fields. The 
experimental data are fitted by exponential curves, in accor- 
dance with equation ([3|. The inset indicates the variation of 
tunnel splitting as a function of the transverse field, the red 
line is a guide for the eye. 



conductance measurements were recorded for a given 
sweep rate. The tunneling probability Pqtm of molecule 
B can be obtained by Pqtm = ^ — Pdt^ where Pdt is the 
probability of a DT [Fig. 1(c)]. Pqtm is plotted in Fig. 6 
where the experimental data were fit with Eq. |3J from 
which a tunnel splitting of A = 1.7 jaK was extracted. 
As a consequence, the mono-exponential behavior clearly 
demonstrates that only the four circled level crossings in 
Fig. 1(d) are avoided level-crossings. This is in agree- 
ment with the work of Vincent et al.^^ but not with the 
measurement done on a single crystal of TbPc2. The lat- 
ter case presents QTM at all intersections of the diagram 
depicted in Fig. 1(d), see the work of Ishikawa et al.^ 
In an assembly of molecules, coupled with weak dipole 
interactions, multi-spin tunnel effects^^^ might be re- 
sponsible for this observation but further investigations 
are needed to better understand this issue. 



VII. LANDAU-ZENER TUNNELING IN A 
SINGLE TbPc2 

An applied transverse field tunes the tunnel splittings 
via the H±{J-^ + J_) term of the Hamiltonian. In order 
to study this effect on the different level crossings, we 
measured the tunneling probability for several constant 
transverse fields and field sweep rates. The symbols in 
Fig. 6 correspond to the experimental points and the con- 



tinuous lines are least-square fits using Eq. (|3|. The data 
agree very well with the Landau- Zener behavior, which 
suggests that no other measurable avoided-level cross- 
ings are induced by the application of a transverse field. 
The tunnel splitting amplitudes were extracted from the 
fits, and then plotted as a function of the transverse field 
(inset of Fig. 6). This behavior cannot be explained by 
using the parameters of Ishikawa et al.^^ Firstly, the lig- 
and field might be different from the mean bulk value 
because the grafted molecules are probably slightly dis- 
torted and their anisotropy modified. ^^ Secondly, the lig- 
and field Hamiltonian does not predict a linear increase 
as observed for our measurements, which were confirmed 
on other molecular devices. This observation needs fur- 
ther experimental and theoretical investigations. In par- 
ticular, we believe that, in the case of single molecules, 
the angular moment conservation has to be taken into 
account. When the latter is not conserved, the Landau- 
Zener equation is not valid. ^^ Nevertheless, our studies 
showed that the hyperfine interaction is a robust feature 
allowing us to read the nuclear spin state regardless the 
deposition and measurement techniques. Moreover, the 
possibility of tuning the tunnel splitting is very conve- 
nient for experiences of coherent nuclear spin manipula- 
tion since the read-out mechanism needs the right value: 
not too small (no read-out) and not too large (DT pos- 
sible, reducing the read-out fidelity). 



VIII. CONCLUSION 

We quantitatively investigated quantum tunneling of 
magnetization at the single-molecule level and confirmed 
that a nanotube-based device is one of the best means 
for the detection of single magnetic moment reversal. 
Furthermore, we demonstrated that the nanotube device 
used for these measurements was able to read-out in a 
single-shot the state of an individual nuclear spin. The 
demonstrated tunability of the tunnel splitting, achieved 
by applying a transverse field, provides the possibility to 
adjust the quantum dynamics of such quantum spintronic 
devices. 
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